Superconductivity and spin-density-waves in multi-band metals 
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We present a detailed description ol two-band quasi-2D metals with s-wave superconducting (SC) 
and antiferromagnetic spin-density wave (SDW) correlations. We present a general approach and 
use it to investigate the influence of the difference between the shapes and the areas of the two Fermi 
surfaces on the phase diagram. In particular, we determine the conditions for the co-existence of 
SC and SDW orders at different temperatures and dopings. We argue that a conventional s-wave 
SC order co-exists with SDW order only at very low T and in a very tiny range of parameters. An 
extended s-wave superconductivity, for which SC gap changes sign between the two bands, co-exists 
with antiferromagnetic SDW over a much wider range of parameters and temperatures, but even 
for this SC order the regions of SDW and SC can still be separated by a first order transition. We 
show that the co-existence range becomes larger if SDW order is incommensurate. We apply our 
results to iron-based pnictide materials, in some of which co-existence of SDW and SC orders has 
been detected. 

PACS numbers: 74.25.Dw,74.25.Ha 



I. INTRODUCTION 

Discovery of new magnetically-active superconduc- 
tors, iron pnictidcs, based on FeAsii 2 - or Fe(Se,S,Te)^ 
has further invigorated the on-going discussions about 
co-existence of different ordered electronic states in 
mctals£~— In itinerant electrons systems, the interactions 
that lead to formation of superconducting (SC) and mag- 
netic spin-density-wave (SDW) orders, "pull" and "push" 
the same particles, and as a result, influence each other. 
In particular, two orders may support each other and lead 
to homogeneous local co-existence of SC and SDW states; 
or one of them may completely suppress the other order, 
resulting in a state with spatially separated regions of 
"pure" SDW or SC orders. The transitions between var- 
ious states may also be either continuous (second order) 
or abrupt (first order). The outcome of this interplay de- 
pends critically on a number of parameters: properties of 
the interactions, such as symmetry of SC pairing, their 
relative strengths, and also on properties the Fermi sur- 
face (FS), such as its shape or the density of electronic 
states. 

In pnictides this parameter space is vast. First, these 
are multi-band materials, with two hole pockets in the 
center, (0,0), and two electron pockets near (±7T, 0) 
and (0, ±7r) points of the unfolded Brillouin zone (BZ) 
(one Fe atom per unit cell). The shapes of quasi two- 
dimensional electron pockets are quite distinct in dif- 
ferent materials, ranging from simple circle-like types in 
LaOFeB^, to cross-like electronic FS in LaOFeAs^ to el- 
lipses in BaFe2As 2 10 ' n and even more complex propeller- 
like structures in (Ba,K)Fe2As2^ (for a descending point 
of view on this see Ref. [Til ). Hole pockets are near- 
circular, but different hole pockets in the same material 
usually have different sizes. 

Second, multiple FSs also create a number of differ- 
ent possibilitic a 13 ' 14 for electron ordering in the form of 



SDW, charge density wave (CDW) states, and various 
superconducting states. The SC states include 1) the 
conventional s ++ -wavc state that has s-wave symmetry 
in the BZ and gaps of the same sign on electron and hole 
FSs; 2) the extended s + ~-state that looks as s-wave from 
a symmetry point of view but has opposite signs of the 
gaps on pockets at (0, 0) and (±tt, 0)r£~— and 3) several 
SC states with the nodes in the SC gap, of both s-wave 
and d-w&ve symmetry^ - — 

As a result of this complex environment, the inter- 
play of magnetic and superconducting orders also shows 
some degree of variations. Most of parent compounds 
of iron pnictides are magnetically ordered. Upon dop- 
ing, magnetism eventually yields to superconductivity, 
but how this transformation occurs varies significantly 
between different Fc-pnictidcs. A first-order transi- 
tion between SC and SDW orders has been reported 
for (La,Sm)Oi_ K F a; FeAs^^ On the other hand, in 
electron-doped Ba(Fei_ a; Co x )2As2 recent nuclear mag- 
netic resonance (NMR)) 26 ' 27 specific heat, susceptibil- 
ity, Hall coefficient,— and neutron scattering experi- 
ments 2 ^ indicate that SDW and SC phases coexist lo- 
cally over some doping range. In the same 122 fam- 
ily, experiments on hole-doped Bai-^K^FcAs^ dis- 
agree with each other and indicate both co-existenc o 30 ' 31 
and incompatibilit y 27 ' 32 ' 33 of two orders. Isovalently 
doped 122 material BaFe2(Asi_ z: P a ;)2 shows the region 
of coexistence 

The goal of the present work is to understand how 
the system evolves from an SDW antiferromagnet to an 
s ++ /s H -wave superconductor, and how this evolution 
depends on the shape of the FS, the strengths of the in- 
teractions, and the structure of the SC order. For this 
we derive and solve a set of coupled non-linear BCS-typc 
equations for SC and SDW order parameters and com- 
pare values of the free energy for possible phases. 

We report several results. First, we find that there is 
much more inclination for co-existence between s" 1 and 
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SDW orders than between the same-sign s ++ -wave state 
and SDW. In the latter case, co-existence is only possi- 
ble at very low T and in a very tiny range of parameters. 
Second, the co-existence region generally grows with in- 
creased strength of SDW coupling relative to supercon- 
ducting interaction. That the co-existence is only pos- 
sible when SDW transition comes first has been noticed 
some time ago ; 37 ' 38 and our results agree with these find- 
ings. Third, when SDW order is commensurate, the co- 
existence is only possible when the following two condi- 
tions are met simultaneously: hole and electron FSs have 
different kp (cross-section areas) and different shapes, 
(e.g., hole pockets are circles and electron pockets are 
ellipses). Even then, SDW and SC orders co-exist only 
in a limited range of parameters and temperatures, see 
Sections [TV] and El Figs. [TO] and [U] below. When SDW 
order is incommensurate, the difference in kp is a suffi- 
cient condition, but again, the two orders co-exist in a 
limited range of parameters/temperatures fFig. 111)). 

We also analyze in some detail the interplay between 
the co-existence and the presence of the Fermi surface 
(i.e. gapless excitations) in the SDW state. The "con- 
ventional" logic states that superconductivity and mag- 
netism compete for the Fermi surface and co-exist if SDW 
order still leaves a modified Fermi surface on which SC 
order can form. We find that the situation is more com- 
plex and the mere presence of absence of a modified Fermi 
surface is not the key reason for co-existence. We show 
that a more important reason is the effective "attraction" 
between SDW and SC order parameters, when the devel- 
opment of one order favors a gradual formation of the 
other order. Specifically, we show that: 

• near the point where the transitions from the nor- 
mal (N) state into SC and SDW states cross, SC can 
develop either via the co-existence phase or via a di- 
rect first order transition between pure SDW and SC 
states, In this range, the SDW order parameter is 
small and SDW state is definitely a metal, Fig. [151 

• at low T, the co-existence phase may develop even 
when SDW state has no Fermi surface (not counting 
bands which do not participate in SDW) . In this situa- 
tion there is no Fermi surface for a conventional devel- 
opment of the SC order, but the system still can lower 
the energy by developing both orders, if there is an 
"attraction" between them. This is the case for 
superconductivity and comparable strength of SDW 
and SC couplings, Figs. [71 [TOTa): 

• the SDW phase at low T can be a metal with rather 
large Fermi surfaces, yet SC order does not develop. 
This is the case when SC order is s ++ , Fig. [lJJ 

The close connection between the co-existence of the 
two states and the symmetry of the SC state has been dis- 
cussed earlier in the context of single-band heavy-fermion 
materials^ This connection gives a possibility to ob- 
tain information about the pure states (e.g., about the 
structure of the SC gap) from experimental investiga- 
tions of the SC - SDW interplay, as it has been recently 




FIG. 1: (Color online) Left: electronic structure of the two 
band model considered in this paper, in the unfolded Brillouin 
zone. The hole FS is in the center, with SC order parameter 
A c , and the electron FSs are at (0, n) and (7r,0), with SC 
order parameter A/. The magnetic order with momentum 
Qo = (0, 7r) hybridize hole and electron FSs separated by Qo, 
but leaves FSs at (±7r, 0) intact. Right: by doping or pres- 
sure one may adjust the size and shape of hole and electron 
bands, and also SDW order parameter can be incommensu- 
rate, with momentum Qo + q. These effects are described by 
FS detuning parameter, <5 kq = [£/(k + q) + £ c (k)]/2. 

suggested. 29 ' 40 

The structure of the paper is as follows. In the next 
section we define the model and derive generic equations 
for the SDW and SC order parameters and an expres- 
sion for the free energy. Then we simplify these formulas 
for the case of a small splitting between hole and elec- 
tron FSs and utilize them in Sees. IIIII through [V] In 
Sec. IIIII we focus on a pure SDW state, with special at- 
tention given to the interplay between ellipticity of the 
FS and the incommensuration of the SDW order. In the 
next two sections we discuss possible co-existence of SDW 
and SC states: in Sec. IIVI we present numerical results 
obtained in a wide range of temperatures and dopings, 
and in Sec. [V] we corroborate this with the analytical 
consideration in the vicinity of the crossing point of SC 
and SDW transitions, and at T = 0. In Scc. lVII wc model 
the case when the splitting between the two FSs is not 
small. We present our conclusions in Sec. IVII1 Some of 
the results reported in this work have been presented in 
shorter publications! 41 ' 42 



II. MODEL AND ANALYTICAL REASONING 
A. General formulation 

Since the basic properties of the SC and magnetic SDW 
interactions and their interplay should not depend on 
the number of bands significantly, we consider a basic 
model of one hole and one electron bands. For pnictides 
this means that we neglect the double degeneracy of hole 
and electron states at the center and the corners of the 
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Brillouin zone, which does not seem to be essential for 
superconductin g 20 ' 22 ! 23 ' 43 " — or magnetic order 42r— 

The basic model is illustrated in Fig. Q] Electronic 
structure contains two families of fermions, near one hole 
and one electron FSs of small and near-equal sizes. Such 
two-band structure yields the experimentally observed 
stripe (tt, 0) or (0, tt) magnetic order which in itinerant 
scenario appears, at least partly, due to nesting between 
one hole and one electron bands, separated by momen- 
tum (tt, 0) or (0, 7r). Other hole and electron bands do not 
participate in the SDW order. We assume that SC also 
primarily resides on the same two FSs, at least close to 
the boundary of the SDW phase. The SC order param- 
eter on the other two bands is not zero, but is smaller. 
Once doping increases and the system moves away from 
SDW boundary, we expect that the magnitudes of the 
SC order parameter on the two electron bands should 
become closer to each other. 

The basic Hamiltonian includes the free fermion part 
"Ho, and the fermion- fermion interactions in supercon- 
ducting and magnetic SDW channels, 



T~L — T~Lo + rlA + Urn ■ 

The free fermion part of the Hamiltonian is 



(2.1) 



Uo = y>(k)4 Q Ck tt + E£/( k ')/k'a/k'a > ( 2 - 2 ) 



k 



k' 



where creation/annihilation c' , c-operators correspond to 
fermions near the central hole pocket (0,0), and /- 
operators describe fermions near the electron pocket at 
Qo = (0, tt) and the fermion dispersions near the pockets 
are 



£ c (k) = [i c 



2m c 



£/(k) 



K y 



2m f x 2m /, 



H (2-3) 



The momenta k are measured from the center of the BZ, 
and k' are deviations from Qo. We assume an inversion 
symmetry, £ c ,/(-k) = £ c ,/(k). 

The pairing interaction consists of many different pair 
scattering terms, but the most important one is the pair 
hopping between the hole and electron pockets ; 13 ' 22 



H A = iXXi,, a ,(k,p) (ctJ_ k pf_ vP ,f pa , (2.4) 

+ /k Q /ik,3 C -P/3' C pa' 



k,p 



For definiteness, we consider SC interaction only in the 
singlet channel, i.e. 

Kifi'a><!*> p) = V & WUW? p>« • (2-5) 

The magnetic interaction between fermions is 
1 x - 



H 



rsdw 

p ' p — k ' k 



^, Q ,(p'p;k,k')x (2.6) 



where we symmetrized the expression with respect to par- 
ticle hopping between (0, 0) — (0, tt) and (0, 0) — (0, — tt) 
pockets for later convenience. We take the interaction 
matrix in a simple form, 

Vffia, (p'p; kk') = Vp^ kk , tr aP ■ <r\ >a . , (2.7) 

with a constant Vp S ,*? kk , = V sdw . 

The evolution of the interaction couplings with energy 
was considered in Ref. [l3l Here we assume that the in- 
teractions for low-energy excitations can be represented 
in terms of fermion couplings to order parameter fields 
in the SC and SDW channels. In the spirit of BCS-type 
approach, we introduce the SC order parameters 

A c (k) Q ^ = E^'p K)k'</-p^/p«'>- 

p 

(2.8a) 

A 7 (k) a/3 = {Ur») afi J^V^ (ia^ic-ve^) , 
p 

(2.8b) 

and the SDW order parameter directed along rh. We 
assume that SDW order parameter has a single ordering 
momentum Q = Qo + q, in which case it is fully specified 
by (m q ) Q/3 = (m q er) Q/3 = m q (mer) Q/3 , where 

p 

P 

(2.8c) 

Since (c pQ ,/ p + q /j) ~ (m q cr) Q( 3, the corresponding elec- 
tronic magnetization, 



m(R) 



p 



4 f 



q/3/ 



- iQR 



is m q cos QR for real m q and is m q cos QR — m q sin QR 
for a complex m q = m q + im'^. In principle, SDW order 
parameter may contain several components with differ- 
ent q, which could give rise to domain-like structures of 
m(R). For recent studies in this direction see Ref. HH. 
We perform the analysis of the co-existence between SC 
order and SDW order with a single q. A more general 
form of the SDW order should not qualitatively change 
the phase diagram for SC and SDW states, however this 
assumption requires further verifications. 

Using the forms of SC and SDW order parameters, we 
write the free and interaction parts in quadratic forms as 
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n ° = \ E [^( k ) °L c ^ + ^(- k ) C -k Q C-ka + £/(k + q) /^ +qa /k +q a + £/(-k - q) /l, 



k-qa/-k-q« 



(2.9) 



= £ [ A c(k)a/3 4a C -k0 + A l( k )a/J C_ kQ C k(3 + A/(k + q) a/3 / k+qQ /l k _ q ^ + A* (k + q) Q(3 /-k-qa/k+q/3 



= ^ E [ m q^/3 /k+qa C k/3 + m q,a,3 C^ ka /- k -q/3 + ™ q af} C kQ / k+q(3 + TO qQ/3 /I k _ qa C_ 



k/3 



(2.10) 

(2.11) 



The Hamiltonian Eq. (|2.1[) can be represented in the matrix form 



ka/3 



V 



&(k) 


Ac(k) ^ 







-A c *(k) itr^ 


-fc(-k) 





m \ (m<r T )^ 


m q (mer)^ 





0(k + q) 


A/(k + q) itr^ 





m q (m<r T ) Q/3 


-A^(k + q) ia% 


-C/(-k-q) , 



with * kQ = (c kQ ,c- kQ ,.^ +qQ ,/- k - qQ ), and * being its 
conjugated column. The two diagonal blocks of the ma- 
trix "H k correspond to a purely SC system with A c and 
Af living on two different bands, and two off-diagonal 
blocks contain SDW field m q that couples fermions be- 
tween the two bands. 

To solve this system of equations for the SC and SDW 
order parameters, Eqs. (|2.8[) . we define the imaginary- 
time Green's function 



G(k,r) Q , = -<T T *(T) ka *(0) k/3 ) = ( G ~ G *1 



which satisfies the Dyson equation, 

G- 1 (k, £n )=ie n -W k , 



(2.13) 
(2.14) 



where e n = nT(2n + 1) are the Matsubara frequencies. 
The system of equations is closed by the self-consistency 
equations for the SC and SDW order parameters in terms 
of this Green's function, 

A c (k) = XX C p T^Tr{(io*)iT + G ff (p, en )} 

P e n 

(2.15) 

A f (k) =J2 V Cp Tj2^{(^ v yf+G cc (P,e n )} (2.16) 

P En 

™ q = -E VSdW jH Tr{(rh* 4 )T 3 G fc (p,ie n )}. 

P En 

(2.17) 

Henceforth we define Pauli matrices in particle-hole 
space, 7^2,3, ?± = (ti ± zf2)/2, and the following ma- 



r 



(2JL2) 



trices in spin- and particle-hole space, 
A = 



o (A i<jy) aP 

(A W)l p 



<T 4 = 



rr \ 

er T ) 

(2.18) 

The expressions above are valid for complex A(k) and 
niq. Below, to simplify formulas, we assume that A's and 
niq are real, i.e., consider only "sinusoidal" , cos QR, vari- 
ations of the SDW order parameter. To lighten the nota- 
tions, we will also drop the momenta arguments (k, k+q) 
in £ c j, A c j and the subscript in m q [still implying this 
dependence as it appears in Eq. f|2, 12[t ] . 

The equations for components of the Green's function 
are obtained from inversion of Eq. (|2,14[) , 



G co ~ m2 Gfo , 

G f } — G fn — m 2 G C Q 



IS 



'jo 



G /c = MG /0 G cc (2.19a) 
(2.19b) 



with definition 
( G-J -M \ . 



ie n - icTz - A c 



-(mcr 4 )f 3 



-(m£r 4 )f3 



ie n - S,fh -A f J 
(2.20) 

To obtain Eq. (|2.19[) we used the fact that the magnetic 
matrix M commutes with purely superconducting parts, 
[M, Geo] = [M, Gfo] = 0, and MM = to 2 . 

The diagonal Green's functions G c o and G/o are the 
same as in a pure superconductor, e.g., 



Gfo(e n ) 



Dfo 



D f o = el+$ + A}, (2.21) 



where for inversion we used the relations 

{f 3 ,A} = 0, [(mo-4)f 3 ,A] =0, A 2 = A 2 , (2.22) 
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which arc also employed to invert 4x4 matrices for mixed 
SC+SDW state. For example, for G cc we have 



and 



1 



G£(e n )-m*GjZ(-e n )/Df ' 
and with the above relations in mind it becomes 



where 



G« 

cc 

q{t 3 ) . 

CC 

£(A) 

cc 

The denominator 



G C c G cc + G cc 3 ^ + G^ c ^ , 
-ie n (D f0 + ni 2 ) 



D 

jcDfO ~ g/TO 2 

A c D f0 - A j 
D 



-73 , 



(2.23) 

(2.24a) 

(2.24b) 
(2.24c) 
(2.24d) 



'kq 



A/Oc + qJ-AcCk) 



A+ _ A f (k + q) + A c (k) 

kq 2 



(2.29) 
(2.30) 



The parameter £k q describes the dispersion and parame- 
ter <5k q describes deviations of the electron and hole FSs 
from perfect nesting, as illustrated in Figs. [1] and [5] 

For the inter-band part of the Green's function we ob- 
tain, 



Gfc{ e n) — MGf G cc 



D 



OD/o + m 2 ) 2 + (ZcD f o - O™ 2 ) 2 + (Ac%o - A f m 2 ) 2 



= MG f0 (e n ) 



G- G l {-e n )D f0 -m 2 G^{e n ) 



D 



D 



fo 



(4 + 4 2 + a 2 )(4 + ^ + a 2 ) 

,2/„2 



= M 



G^(-e n )G^(-e n )-m 2 
D . 



2ml4 - CcC/ - A c A f ) 



(e 2 n + E 2 + )(e 2 n + E 2 _) 



(2.25) 



(2.31) 

where for self-consistency Eq. (|2 . 1 T[) we need only the 
part proportional purely to M-matrix, Eq. (|2.20l) . 



gives the energies of new excitations in the system, c.f. 
Ref. [29|. We obtained (explicitly showing k and q here): 



MM) _ , ^ -sl + j fZc + A f A c 



El = CL+^ q + m 2 + (A kq ) 2 + (A+) 



(2.26) 



D 



(2.32) 



± 2 x /m 2 [(A+ q ) 2 + <5 2 q ] + (A" q A+ q + £ kq «i kq ) 



kq kq 



with 



£kq 



^kq 



Q(k + q)-g c (k) 
2 

Q(k + q)+^(k) 



(2.27) 
(2.28) 



Expressions for G// and G c / are obtained from Eq. (|2.24[) 
and Eq. f|2.32[) by swapping indices, cf>/. 

We substitute the above expressions for G c f\p), 
Gff (p + q) and GrM' (p) into the self-consistency equa- 
tions Eqs. ((2~T5l - PTT71) and arrive at 



A / (k)=T^^(-2V t 



k p) 



A c (p)[e 2 + (p + q) + A 2 (p + q)] - A/(p + q) 

D 

~D 



A c(k) = rS£(-»g^ A/(p + q)[4 + eAp) D + Ac(p)] - Ac(p) 

e„ p 

to _ t ^^ t , s a u 4 + ^ 2 -Q(p + q)6(p) - A / (p + q)A e ( P ) ^ 



(2.33) 



(2.34) 



(2.35) 



To calculate the relative stability of different states, one also needs to evaluate the free energy. We follow the 
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Luttingcr-Ward 5 ^ and Dc Dominicis-Martin 5 ^ method, 
and consider the functional^ 

F = -lsp {EG + In[-(ie„ - £) + E]} + $[G] , (2.36) 

which, if minimized with respect to G, gives self- 
consistency equations, E[G] = 2<5$[G]/<5G; and, if min- 
imized with respect to E, gives the Dyson equation, 
Eq. (|2.14|) . Here Sp is the trace over two fermion bands, 
spin, particle-hole matrix structure, and the sum over 
Matsubara energies and the integral over momenta, and 
E is the mean held SC and SDW order parameter matrix, 

£ -(*£)- 

The functional $[G] producing the self-consistency equa- 
tions is a quadratic function of G. Using the self- 
consistency equations one can explicitly verify that at 
weak-coupling it can be written as $[G] = |Sp{EG}. 
To deal with the logarithm in Eq. (|2.36|) one introduces 
a continuous variable A instead of e„, differentiates the 
logarithmic term with respect to A to obtain the Green's 

function G(A) = (iA — | — E^ , and then integrates 

back to get the difference between a condensed state and 
the normal state for fixed external parameters, such as 
temperature or field, 

AF(A c j, to) = -^Sp | i£G - Jd\[iG(X) - iG N (X)} J 

(2.38) 

where Gm is the Green's function in the normal state 
without either SC or SDW order, and we used the fact 
that in the normal state E = 0. Substituting into 
([235)1 the Green's functions Eqs. (j2~2"4"[) . (g32J) , the self- 
energy Eq. (|2.37[) , and using the self-consistency equa- 
tions Eqs. (|2.33[) - (|2.35[) to eliminate the high-energy cut- 
offs in order to regularize the e„-summation and k- 
integration, one obtains the most general free energy 
functional for given A c j and to. 

B. Limit of small Fermi surface splitting 

In principle, equations for full Green's functions 
Eq. (|2.24[) . (|2.32l) . the self-consistency equations 
Eqs. ([2~53 > ([2~31) ]l and the free energy Eq. ([2~55)l . com- 
pletely describe the system in a very general case. How- 
ever, to proceed further with the analytics one can reduce 
the number of summations which is also desirable from 
a numerical standpoint. 

The typical approximation is to linearize the dispersion 
near the FS and integrate out the momenta in the direc- 
tion normal to the FS over £kq- In the case, when the 
two FSs are reasonably close to each other (when shifted 



by (0,7r)), and electron and hole dispersions are similar, 
the values of FS mismatch (5k q are weakly momentum 
dependent, and can be taken at positions where £k q = 0. 

The consequence of this approximation, which we dis- 
cuss in some detail in Appendix [XI is that (5k q depends 
only on the angle in k-space, but not on £kq and hence 
one can integrate along a particular direction k over £kq, 
keeping J^. fixed. 

Within this approximation the DOS for both FSs are 
the same, and the magnitudes of A c and Af are equal 
(the angular dependence of SC gaps is still determined 
by that of the SC interactions). There are, indeed, also 
higher order terms, which we neglected in the last lines 
of Eq. (jA2[) . These terms make hole and electron DOS 
different from each other, what in turn makes |A C | and 
| Af | non-equal, but these terms are small in <5k//-t c ,/ an d 
only account for sub-leading terms in the free energy, jj, c j 
are Fermi energies of electron and hole bands, Eq. (|2.3[) . 

This approximation comes at certain price. When two 
FSs are of very different shapes, approximating them as 
small deviations from a single line in k-space everywhere 
is incorrect. This is shown for example in Fig.[2jd), where 
the two FSs are quite different away from the crossing 
points. However in this case one realizes that if at some 
k-point the two bands arc far apart, the effect of the 
SDW is very small, and we can approximate those FS 
parts as participating in SC pairing only, with little or no 
competition from the SDW interaction. This can be seen 
from Eq. (|2.19p for the Green's function. For example, 
for electrons near the FS of the c-band, £ c — > 0, £/ is large 
and G~} ~ G~q + 0(to. 2 /£/), and the corrections due to 
m can be neglected when we go along c-FS away from the 
region where £ c ~ w 0. We will return to this issue in 
section fVTl to show that the results are qualitatively the 
same whether we consider large or small splitting of the 
FSs. 

For small splitting between hole and electron Fermi 
surfaces, we perform ^-integration analytically. For this 
we approximate V£ c by an isotropic V sc , i.e., take angle- 
independent SC gap. The sign of V sc can be arbitrary, 
and we consider separately the two cases: 

a) V sc > 0: results in the s + ~ state, with gaps of 
opposite signs for electrons and holes, 

A/ = -A c = A or A + = 0, A_ = A; 

b) V sc < 0: s ++ state, with the same gaps on two FSs, 

A f = A c = A or A + = A, A_ = . 

In both cases A + A_ = and the denominator of the 
Green's function can be written as, 

D = (e 2 l + E 2 + )(e 2 n +E 2 _) = (^ q + E 2 + )(£ 2 q + E 2 _) , (2.39) 

where 

El = el+A 2 +m 2 -5 2 ±2\ m 2 A 2 ^^- - 5 2 (e 2 + A 2 ) 

it n kqy 2 kq vn ' 

(2.40) 
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with s = +1 (s 
state. Closing the integration contours over £k q in the 



-1) corresponding to 



self-consistency equations and in the free energy over the 
upper half-plane and counting poles at +iT,± wc obtain 




AF(A,m) A 2 , T m 2 , T v-^ / 1 , . A 2 m 2 

K —± — - = — n — + — In 2ttT > ( -(£+ + E_) - e„ 

4iV P 2 T, 2 T, ^ \2 V + ; 2e„ 2e„ 



En>0 

r 



In- 



sdiv 



1.13A 



In 



1.13A 

T s 



(2.41) 

(2.42) 
(2.43) 



where angle brackets denote remaining momentum av- 
eraging over directions on the FS. Np is the density 
of states at the FS per spin, and v sc = 2NpV sc and 
y sdw _ jy^ysdw are dimensionlcss couplings in the 
SC and SDW channels^ Taken alone, v sc leads to a SC 
state with critical temperature T c , which is independent 
of <5^ q as one can check by setting m = in Eq. (|2.4ip . 

while v w leads to an SDW state with transition tem- 
perature T s which does depend on 8^ . We define T s as 
the SDW transition temperature at perfect nesting, when 

The relative sign between SC and SDW orders, as 
given by terms -sAm 2 in Eq. (|2.41[) and — smA 2 in 
Eq. (|2.42p . is positive for s H state resulting in effec- 
tive "attraction" of the two orders, and negative for s ++ 
state implying that the formation of one order resists 
the appearance of the otheri 29 i 41 The actual co-existence 
of the two orders, however, is a more subtle effect and 
needs to be determined from the exact solution of these 
equations and the analysis of the free energy. The dif- 
ference in excitations energies Eq. (|2.26p and Eq. (|2.40[) 
between s ++ and s" 1 states is also consistent with previ- 
ous studies of d- and p- wave superconductivity in heavy- 
fermion metals^S that concluded the SC states with sym- 
metries PTq = — 1 (e.g. s" 1 , d), where P is parity 
[PA(p) = A(— p)] and Tq is the shift by the nesting 
vector [T<gA(p) = A(p + Q)], are more likely to form 
co-existence with SDW than those with PTq = +1 (e.g. 

8++,p). 

We will also analyze the quasiparticle density of states 
(DOS), which is given by the integrals over £kq of the di- 
agonal components of the Green's function. For example 
for c-fcrmions 



9c{e n , k) 



-is Tl 



G (1) 



(2.44) 




which for pure SDW state reduces to 



9c {e 



».*) = -;£ 



2 ± Jm 2 - {ie n ±S Uci ) 2 



(2.45) 



and actual DOS is obtained by analytic continuation 
jV(e,k) 



N F 



-3mg(ie n -+ e + iO+,k). (2.46) 



III. PURE SDW STATE 

In pnictides, parent materials usually have only mag- 
netic order below a transition temperature T s . Super- 
conductivity appears at a finite doping, when the SDW 
transition is suppressed. Keeping this in mind, we con- 
sider first a purely SDW state, and analyze how it is 
modified when FSs are deformed by addition or removal 
of electronic carriers, and whether modified FSs are still 
present in the SDW phase. To remind, we denote by T s 
the SDW-N transition temperature at perfect nesting, 
which effectively gives the scale of SDW interaction in 
the system. The true instability temperature, which we 
denote explicitly by T s (6? ), is a function of ellipticity, 
doping, and incommensurability. 

We begin by presenting explicit formulas for the exci- 
tation spectrum, the SDW order parameter, and the free 
energy. For A = 0, 



Y? ± given by (pOO]) is 



S 2 



(en ± % q 



(3.1) 



and the excitation spectrum consists of four branches 
with energies ±E±(A = 0), where 



E±(A = 0) 



C 2 q + m 2 ±8 



kq 



(3.2) 
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In JO and (0L2 



a) 



2 

£/(k + q)+£ c (k) 0(k^ + q) 



k q 2 2 

= ^ k (k^-4 + q). 



(3.3) 



We remind that 5/ describes the mismatch between the 
shapes of the electron and hole bands and determines 
their nesting properties in k-direction. 

Equation (|2.42[) for the SDW order parameter to sim- 
plifies to 



4- = 2ttT V Re 

saw / -j 



l 



£„<A \ ( £ n + iS i a ) 2 + m2 



-, (3.4) 



and the cut-off A can be eliminated in favor of T s , 

T / 1 1 \ 

In — = 2ttT V Re ( — -. — r ) , 

Ts £U \ V /(^+*^ q ) 2 +m 2 M/ 

(3.5) 

where the summation over e„ now extends to infinity. 
Second-order transition temperature T — T S (S^) is ob- 
tained by setting m = 0: 



ln^=27rT J2 Rc 



e„>0 \ £ «+%q £ - 



(3.6) 



The free energy, Eq. (|2.43|) . becomes 

AF(to) _ to 2 T 
4iV F ~ "2" n T7 



£„>0 

to 2 , 1.13A 
m — — — 

2 T, 



-2^T ^ Re / J(e n + z% q ) 2 + m 



m 
2^ 



0<e„<A 



2^T ^ (Re U(e n + iS^) 2 +mA-e n ) (3.7) 



Below we consider several special cases for 5^ (see 
Fig. [2): 

• two co-axial circles, q = 0: — k£ = k.(kp — kp): 



(3.8) 



For a fixed 5q, circular hole and electron FSs survive in 
the SDW phase when m < 5$ (Fig. Wp), but come closer 
to each other as to increases and merge at to — 8q. At 
larger to all excitations are gapped (Fig. [3h). 

• FS of different shapes, e.g., one circle and one el- 
lipse, co-centered: q = with k^ — k£ = \t(kp — k F + 
Ak cos2</>): 



% q = 5q + 82 cos 20 , S 2 = -v F Ak (3.9) 



■l: f . -i; 

v. 

'•>•'••. . 




b) 




d) 



FIG. 2: (Color online) The appearance of gapless excitations 
in the presence of SDW order. The dotted lines indicate FSs 
for electrons, — 0, and holes, £ c = 0. The dashed curve is 
an "effective" FS, £k q = 0. a) when q = and m is large com- 
pared to FS mismatch, m > 5^ , all excitations are gapped; 
b) when m is small, gapless excitations are preserved along 
the two modified FSs at £k,o = i(<^ — m 2 ) 1 ^ 2 Cl^ft 1 > m i n 
the shaded region). Such gapless state, however, only exists 
at high temperatures, while at low T it is pre-emptied by a 
first order transition to the normal stated; c) to prevent the 
first order transition, magnetic order is formed at an incom- 
mensurate vector Q = Qo + q. This improves electron-hole 
nesting on some part of the FS, but allows for gapless excita- 
tions at the opposite side; d) when the two FSs are of different 
shapes, the nested parts become gapped due to SDW order, 
and on the rest of the FSs the excitations are little affected 
by SDW order. The density of states for these cases is shown 
in Fig. El 



In this case, at small enough to, the FS has a form of 
two hole and two electron pockets. As m gets larger, the 
pockets shrink and eventually disappear. 

• two circles of different radii, centers shifted by q: 



^kq = S U = S ° 



^v F q = 5 Q + ^v F qcos((j) - <f> ) (3.10) 



where (f> and 4>o are the directions of v F and q. In this 
case, when m increases, gapless excitations survive along 
a pocket in one region of the k, while excitations with 
— k become gapped ( Fig. [2b). At large enough to, mod- 
ified FS disappears and excitations with all momenta k 
become gapped. This scenario refers to the case when 
the magnetic ordering occurs at a vector, different from 
the nesting vector Qo, producing incommensurate SDW 
state. It may occur because the electronic system has an 
option to choose q 7^ if it minimizes the energy, or be- 
cause the SDW interaction is peaked at a fixed Q 7^ Qo 
for some reason. Note that Eq. (|3.5|) for SDW order 
is a magnetic analog of Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) stated in a paramagnetically limited supercon- 
ductor. An incommensurate SDW state at finite dop- 
ings has been studied in application to chromium and its 
alloys^— and, more recently, to pnictide materialsi 42 i 60 



9 



2\ 



1 


i 1 i 1 




_ \ = 0.13 


\ 


0.00 

0.05 

— 0.10 
0.15 






25, 























O.I 



0.2 0.3 

EllTZT 



0.4 



0.5 



FIG. 3: (Color online) The FS averaged DOS for pure SDW 
state and FS mismatch 8^ — So + 82 cos 2cj>. We use di- 
mensionless variables denoted by bars, 80 = So/2irT s = 0.13, 
zero-temperature SDW gap mo = m(T — 0)/2tyT s = 0.28 
and vary 82 = 82/2nT 3 . For 82 = 0, N(e) vanishes below 
e = m — 80 and has two sharp BCS peaks at e = m ± So- At 
finite ($2, each of the two peaks splits into a "band" bounded 
by two weaker non-analyticities separated by 282 . The gap in 
the DOS behaves asm- 80 — 82 and closes when 80 + 82 > m, 
metallic states forms. The DOS remains the same if we re- 
place the ellipticity parameter 82 by the incommensurability 
parameter q. 



In general, all three terms are present, and 

d kq = S o + S 2 cos 20 + ^v F qcos{(/) - O ) ■ (3.11) 

In the figures we use dimensionless parameters, that are 
denoted by a bar. For isotropic and anisotropic FS dis- 
tortions, 



r t>0,2 _ v F q 

00,2 = TTT^- j 1 



2nT s ' 4ttT s 1 

and similarly for other energy variables, 
m e T A 



2ttT, 



2ttT, 



2ttT c 



(3-12) 



(3.13) 



We use different notations for prefactors of cos(0 — 0o) 
and cos 20 terms to emphasize that they have different 
origin: S 2 is an "input" parameter defined by the elliptic 
form of the electron FS due to the electronic band struc- 
ture, while q is adjustable parameter that minimizes the 
free energy of the system. If the minimum of the free en- 
ergy corresponds to q = 0, SDW order is commensurate, 
otherwise SDW order is incommensurate. 

In Fig. we show the DOS N(e) for the fixed S = 
0.13 and m = 0.28, and different 8 2 . For S 2 = 0, N(e) 
vanishes below e = m — 5q and has two BCS-like peaks 
at e = m±5a. At finite S 2 , each of the two peaks spreads 
into a region of width 28 2 bounded by two weaker non- 
analyticities. The gap in the DOS behaves as m — Sq — 8 2 
and closes when Sq + 5 2 become larger than m. The 
DOS and all other results remain the same if we replace 
the ellipticity parameter 5 2 by the incommensurability 



parameter q because the angular integral in Eq. (|2.45[) or 
Eq. (|3.6p over momentum directions on the FS coincides 
for cos(0 — 4>o) and cos 20 terms in 5^ , if considered 
separately. The DOS and ^(5^) change, however, when 
both 5 2 and q are present simultaneously. 

Below we discuss the phase diagram for the pure SDW 
state to the extend that we will need to analyze potential 
co-existence between SDW and SC states, which is the 
subject of this paper. 

It is instructive to consider separately the case when 
SDW order is set to remain commensurate for all Sq, 2 
(i.e., q = 0), and the case when the system can choose 
q. In our model, the first case is artificial and just sets 
the stage to study the actual situation when the value of 
q is obtained by minimizing the free energy. However, a 
commensurate magnetic order may be stabilized in the 
SDW state, if the interaction V sdw is by itself sharply 
peaked at the commensurate momentum Qo. 

The results for the case q = are presented in Fig. [4] 
In panel (a) we present the results for the transition tem- 
perature T s (5o,5 2 ) for several values of 5 2 . All curves 
show that the transition is second-order at high T and 
first-order at small T. The first-order transition lines 
(dotted lines in Fig. SJa)) were obtained by solving nu- 
merically the nonlinear equation for m, substituting the 
result into the free energy (|3.7[) and finding a location 
where AF(to) = 0. 

To verify that the transition becomes first order at low 
T, we expanded the free energy in powers of m as 



AF(ra) = a n 



Bm 4 



(3.14) 



and checked the sign of the B term. The coefficients a„ 
and B are determined from Eq. (|3.7I) . 



a m = ~ (in J- -2nT ^ Re 

\ 3 e n >0 



1 1 



i6 U 



4 e ^cT~ \( £ » + %) 3 



(3.15) 



where 5^ = Sq + S 2 cos 20. Solid lines in Fig. 0|a) corre- 
spond to a, n = 0. The N-SDW transition is second order 
and occurs when a, n = if B > 0, but becomes first 
order and occurs before a m becomes negative if B < 0. 
We indeed found that for all fixed 5 2 , for which SDW- 
N transition is possible, B changes sign along the line 
a m = and becomes negative at small T. For S 2 = 0, 
this occurs at T* = 0.56 T s and 8% = 0.17. 

We point out the following counter-intuitive feature in 
Fig. |4|a). Increase in 6 2 reduces the transition temper- 
ature at 5o = 0, and at the same time makes the curve 
flatter allowing for a larger SDW region along So. The 
transition line becomes completely flat at a critical value 
5 2c = 0.28073(2ttT s ) (see below) when T S (S , S 2c ) = +0. 
At this point, it spans the interval 5q £ [0, S 2c \. The exis- 
tence of the SDW ordered state at = S 2c over a finite 
range of So despite that the transition temperature is +0 
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FIG. 4: (Color online) The SDW-N transition for commen- 
surate SDW order. The parameters So and 82 describe the 
difference between the area of hole and electron pockets and 
the ellipticity of the electron pocket, respectively. Here and in 
all subsequent figures dotted lines mark first-order transitions, 
solid and dashed lines mark second-order transitions. Panel 
(a): variation of the transition temperature with So for fixed 
&2- The transition is second order at small So but becomes 
first order at larger So. At S2 = 0, the transition becomes first 
order at T* « 0.56T S . Panel (b): variation of the transition 
temperature with 82 for fixed So- T s (So,S2) monotonically 
decreases with increasing S2 and vanishes at the same value 
82 ~ 0.28073, independent of S . 
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FIG. 5: (Color online) Same as in Fig. [4] but when the sys- 
tem is allowed to choose between commensurate and incom- 
mensurate SDW orders. Solid lines are second-order tran- 
sition lines into a state with a commensurate SDW order, 
dashed lines are second order transition lines into an SDW 
state with an incommensurate SDW order (the magnetic ana- 
log of FFLO state). For all £2 > 0, incommensuration occurs 
before the commensurate transition becomes first order (the 
onsets of incommensuration and first-order transition coin- 
cide for 82 = 0). Observe that incommensuration develops at 
progressively smaller So as 82 increases and T 3 (So) decreases, 
but the range of So over which incommensurate SDW order 
exists actually increases with increasing 82- 



is a highly non-trivial effect which deserves a separate 
discussion^ 

In Fig. SJb) we show the transition temperature at 
fixed 5q, as a function of the ellipticity parameter 82 ■ As 
expected, Tg^a) monotonically decreases with increasing 
ellipticity of the electron band. The SDW order exists 
up to 62c, at which T s (d2 C ) = +0. The value of 62c is 
independent of So and can be obtained by taking the 
limit T — > in (|3.4p with m = and re-writing this 
equation as 



1 f I 1 

—r~ = Re ds 



e + iSfc 



Re In 



2A 



i6 + - 81 



(3.16) 

The interaction can be eliminated in favor of zero- 
tcmpcrature gap mo at 5q = 82 = 



1 



,sdw 



de- 



1 



In 



2A 



^/e 2 + ml mo ' 



(3-17) 



where from Eq. (|3.5p we obtain, at So = 82 = 0: 
2wT s 



m 



0.28073 x (27rT s ) 



(3.18) 



and 7b w 0.57722 is Euler's constant. At finite 80 and 
8 2 , the value of m at T = remains equal to m a as long 
as 80 + 82 < m - The combination of Eqs. (|3.16[) and 
(t3~T7)) gives 8 2c =m = 0.28073 x (2ttT s ), provided that 



Sq < S2C (there exists another solution S\ = 2Somo — m\ 
at 82 < So and So > mo/2, but it corresponds to an 
unstable state). A similar result has been obtained in 
the studies of FFLO transitions^ 

The form of T S (S2) near 82c can be obtained analyti- 
cally by re-writing the condition a m = in (|3.15p as 



In 



T 



2nTj2 



£n>0 



82 cos 2 2(j> 



Sn(el + Si cos 2 2(j>) 



0, (3.19) 



integrating explicitly over <fr, re-expressing 1 / \/sf t 
as (2/n) / °° dxj (x 2 + e 2 + <5f ), and performing the sum- 
mation over e n before the integration over x. Carrying 
out this procedure, wc obtain 



T s (5 2 ) 



'2c 



ln(l-5 2 /(5 2c )| 



(3.20) 



We see that T s very rapidly increases at deviations from 
S 2c . For S 2 = 0.99745 2c (S 2 = 0.28), Eq. (I3~2U)) yields 
T S {S2) ~ 0.3T S , in good agreement with Fig. B|a). 

Also, one can easily show that at T = fcrmionic 
excitations in the SDW state are all gapped when mo > 
80 + 82- When S0 + S2 > mo, the SDW state possess Fermi 
surfaces and gapless fermionic excitations. 

We next consider the case when the system is free 
to choose between commensurate and incommensurate 
SDW orders and may develop incommensurate order to 
lower the free energy. In Fig. [S] wc show the transition 
temperature T s (5o) for fixed 82 ■ We found that, for all 
82 , first order transition is overshadowed by a transition 
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into an incommensurate SDW state. For 62 = 0, incom- 
mensuration develops exactly where B changes sign, and 
the transition into incommensurate SDW state remains 
second order for all Sq. For 82 > 0, incommensuration 
develops before B changes sign, and the transition into 
incommensurate SDW state remains second order over 
some range of 80 but eventually becomes first order at 
large 5q and low T . The full phase diagram also contains 
a transition line (not shown in Fig. [5]) separating already 
developed commensurate and incommensurate SDW or- 
ders. 

To analyze the interplay between the appearance of 
incommensurate SDW order and the sign change of B, 
we again expand the free energy in powers of m but 
now allow incommensuration parameter 8\ to be non- 
zero, i.e., replace in the coefficients in Eq. (|3.14p , 8^ = 
So + 62 cos(2</>) with 8^ = 8^ + q cos(4> — (fio). In general, 
for small q, 



(<*k, q ) = MSi) + a2(S t )q 2 +a 4 (S t )q 4 +0(q (i ), (3.21) 



with ao^k) gi ven by (|3.15p . When 014 and B are posi- 
tive, the N-SDW transition is second order, and is into 
a commensurate SDW state when a 2 > and into an 
incommensurate SDW state when 0:2 changes sign and 
becomes negative. If B changes sign while ct2 is still pos- 
itive, the SDW-N transition becomes first order before 
incommensuration develops. 

To understand the phase diagram, it is sufficient to 
consider small 82 ■ Expanding all coefficients in powers of 
82 we obtain 



012(61) 

a 4 (<y 

B 



a ,o + a ,2 <5f + 0(<£) , (3.22a) 
a 2 ,o + a a ,i cos 2</> 82 + 0(5l) , (3.22b) 

o(s 2 2 ), 



"4,0 
1 



00,2+00*1). 



(3.22c) 



where 



"2,0 




a 2 ,i = -2ttT 



e„>0 



f>4 



16 



2^E 



108 



1J2 
b n 



(4 + ^o 2 ) 5 



(3.23d) 



We see from Eqs. ([3321) that for 8 2 = 0, B and a 2 (6 is ) 
change sign simultaneously, at the point where 0(2.0 = 
a o,2 = 0. However, when 82 7^ 0, 0.2 (6^) changes sign be- 
fore B becomes negative because 0:2(6^) contains a term 
linear in 82, whose prefactor can be made negative by ad- 
justing <j)Q. This explains why in Fig.[5]incommensuration 



begins while B is still positive. Also, we verified that near 
the onset points for incommensuration, a±(8^) > 0, i.e., 
in this range the transition into incommensurate SDW is 
second order. At larger 80, the incommensurate transi- 
tion eventually becomes first order. 



IV. SDW+SC STATE, NUMERICAL ANALYSIS 



In the next two sections we look at potential co- 



existence of SDW and the 



or s ++ states, when the 



system is doped and the SDW state is suppressed. The 
superconducting T c is doping independent, so at some 
doping SDW and SC transition temperatures cross. Near 
this point, the two orders either support or suppress each 
other and either co-exist or arc separated by a first-order 
transition. 

In this section we present numerical results in the ex- 
tended range of temperatures and dopings, in the next 
section we corroborate them with analytical considera- 
tion in the vicinity of the crossing point, when both order 
parameters arc small, and at T = 0. 



A. Coexistence with s state 

We look first at the s H state. In this case the sys- 
tem of coupled self-consistency equations for A and m 
is obtained from Eqs. (|2~41^ - (|2~i3"f by taking T,% = 

\2 



and E n = y/el + A 2 



T 



In- =2^^ Re 

C e„>0 



(E n + i6^)/E n 



(£»+% q ) 2 



(4.1a) 




FIG. 6: (Color online) The phase diagram of SDW and SC 
s H states when only a commensurate SDW order is allowed 
(q = 0). We set T 3 /T c — 3 and varied either the relative 
radius of circular hole and electron pockets (a) or the form of 
one of the pockets (b). The pure SC s" 1 and SDW states are 
separated by first order transition, and there is no co-existence 
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FIG. 7: (Color online) Appearance of co-existence when both 
Sq and 82 are finite. We set T s /T c = 2 and q = 0. Panels 
(al)-(a3) - phase diagrams in variables T, So at fixed 82, pan- 
els (bl)-(b3) - phase diagrams in variables T, S2 at fixed So- 
Panels (al), (bl) - there appears a region near T = 0, where 
SDW and SC s H orders co-exist. Panels (a2),(b2) - the co- 
existence region broadens and reaches T — T c . Panels (a3), 
(b3) - the transition at low T becomes first order between 
pure SDW and SC states, but narrow co-existence region is 
still present near T c . A complimented zero-temperature phase 
diagram is presented in Fig llOl 



In- 



T 
T~ s 



2irT Re 



e n >0 \ J {E n T ^ kq ; 

(4.1b) 

We remind that T c is the transition temperature for the 
pure SC state, and T s is the transition temperature for 
the pure SDW state at 5 kq = 0. 

These equations are solved numerically to find all pos- 
sible states (A, m) and their energies evaluated using 
Eq. (|2.43p . The main results for this part are presented 
in Figs. IMH1 



1. Commensurate SDW state 

Figure [6] shows the results for the case when SDW or- 
der is set to be commensurate (i.e., q = 0) and the FSs 
are either co-axial circles (panel a), or of different shapes 
with equal kp (panel b). In the first case, 82 = and 
So ^ 0, in the second case So = and S2 5= 0. We see 
that in both cases pure SDW and SC states are separated 
by a first-order transition. We verified that in both cases 
fcrmionic excitations in the SDW state are fully gapped 
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FIG. 8: (Color online) (a) Same as in Fig. 0F>) but for 
T s /T c = 5. (b,c) SDW and SC gaps, in units fh and 
A = A/2ttT c , and the free energy as functions of So along 
the line T/T c = 0.2 (b), and as functions of temperature for 
So = 0.17 (c). 



at T = and thus there are no Fermi surfaces. From this 
perspective, the results presented in Figure [5] are consis- 
tent with the idea that co-existence requires the presence 
of the Fermi surfaces in the SDW state. However, we will 
see next that the situation in the cases when both So and 
S2 are non-zero is more complex. 

This is demonstrated in Fig. [7] which shows the phase 
diagram for T s /T c = 2 as a function of So for a set of 
fixed S2 (panels (al)-(a3)), and as a function of S2 for 
a set of fixed So (panels (bl)-(b3)). For all cases, pure 
SDW state is fully gapped at T = 0, so naively one should 
not expect a co-existence state. However, as is evident 
from the figure, the phase diagram does involve the co- 
existence phase, which can be either at low T (including 
T = 0), or near T = T c , depending on the parameters. In 
particular, as 82 in panels (a) or So in panels (b) increase, 
the co-existence state first appears at low T, while at 
higher T the pure SDW and SC states are still separated 
by first-order transition (panels (al) and (bl)). Then 
the co-existence region grows, and extends up to T = T c 
(panels (a2) and (b2)). At even larger 82 or So, SDW 
and SC states are separated by the first-order transition 
at low T, but the co-existence phase still survives near 
T c . 

In Fig. [8] we show the phase diagram for 82 = 0.2 
and T s /T c = 5 together with the plots of SDW and SC 
order parameters and the free energy. We see the same 
behavior as in Fig. [7] (a2) - there is a co-existence phase 
for all T up to T c . In Fig. [9] we show the changes in the 
quasiparticlc DOS at low T = 0.1T C as the system evolves 
from the SDW state to the SC state via the co-existence 
region. 

Finally, in Fig. llOl we show the zero-temperature phase 
diagram in variables 80 and 82 for T s /T c = 2 and T s /T c = 
5, together with the locus of points where T s (Sq, 82) = T c . 
The phase diagram was obtained by numerically solving 
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FIG. 9: (Color online) FS averaged DOS as a function of en- 
ergy at T/T c = 0.1 for different values of 5q (i.e., different dop- 
ings). We set 82 = 0.2 and T 3 /T c = 5. The characteristic val- 
ues of the SDW order parameter for this range of parameters 
is m m 0.28 (zero-T limit). Upper panel: DOS for small So, 
when the system remains in the pure SDW state. This figure 
is similar to Fig. [3] The DOS vanishes below e = m — 80 — 82 
and has In — non-analytic behavior at e = m ± 60 + 62, and 
sudden drops at e = m±<5o — 82- Lower panel: DOS for larger 
80, when SDW and SC orders co-exist. Sharp peaks at small 
e are due to opening of the superconducting gap A. Once 
SDW order disappears at 80 ss 0.21, the DOS acquires BCS 
form with the maximal gap A m(l;r /27rT c « 0.28. 



Eqs. (|4.1[) and evaluating the free energy at T/T c = 0.02. 
This phase diagram corroborates the results of Fig.[7Jand 
[5]- the zero T behavior in panels (a) and (b) in Fig. UJ) 
is obtained by taking either horizontal or vertical cuts 
in Fig. [TUJ In particular, we see from Fig. QJJ1 that the 
transformation between panels (a2) and (a3) of Fig. [7] is 
such that the co-existence region at T = first moves to 
the left, shrinks, and disappears at 82 ~ 0.24. Similarly, 
in panels (b), the co-existence range shrinks to a point 
at So f=a 0.16, and at larger So the transition between 
SDW and SC phases at T — becomes first order. In 
the next section we present the results of complimentary 
analytical studies of the phase diagram at T = and 
near T c . These results arc in full agreement with the 
numerical analysis in this section. 

Observe that for T s /T c = 5 the left boundary of 
the co-existence region is located very close to the line 
S + 6 2 = m (dashed line in Fig. [TUl S + S 2 = 0.28073), 
at which gapless excitations and Fermi surfaces appear 
in the SDW state. For this T s /T c , the co-existence region 
at T = virtually coincides with the region where SDW 
state has a Fermi surface. However, for smaller T s /T c — 2 
(Fig. [7J left panel in Fig. [TU|) , co-existence clearly occurs 
already in the parameter range where SDW excitations 
are all gapped. The co-existence for T s /T c = 2 is there- 
fore not the result of the "competition for the Fermi sur- 
face", but rather the consequence of the fact that the 
system can gain in energy by reducing the SDW order 
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FIG. 10: (Color online) The zero temperature phase diagram: 
SDW, SC states and their co-existence region for various 82 
and 8 and two different T s /T c = 2 (left) and T s /T c = 5 
(right). We only allow the system to develop a commensurate 
SDW order (q — 0). The dashed line denotes first appearance 
of gapless excitations in the pure SDW state [m = So + 82, 
Eq. (|2.45[) ]. The co-existence region at T = (shaded area) 
extends down to So = 0, but is not present at small 82- The 
width of the co-existence region increases with the relative 
strength of SDW interaction, as determined by ratio T s /T c . 
The squares mark the location of the crossing point between 
T a (So,S2) and T c . Open squares indicate that the SDW-SC 
transition near T c is first order, while filled squares signal the 
presence of the co-existence phase near T c . Note that the 
regions where co-existence phase is present at T = and near 
T c are not identical. 

parameter (still keeping all fermionic excitations gapped) 
and creating a non-zero SC order parameter. The gain 
of energy in this situation can best be interpreted as the 
consequence of the attraction between the two orders. 



2. Commensurate vs. incommensurate SDW state 

One of the results of our consideration so far is that, 
if we keep an SDW order commensurate, a finite region 
of SDW + SC phase appears only when both So and 
S2 are non-zero. If we allow the system to choose the 
ordering momentum of the SDW state, the co-existence 
region widens and appears even if we set 82 = 0. We 
illustrate this in Fig. [Til where we plot the phase diagram 
at 82 = for two different values of T s /T c . In agreement 
with Fig. El at T < T*, the system chooses an SDW 
state with a non-zero q. We see that, in this situation, 
there appears a region where SC state co-exist with an 
incommensurate SDW state.— The co-existence region 
widens up when the ratio T s /T c increases, and for large 
enough T s /T c extends down to T = 0. In Fig. [12] we set 
S2 to be non-zero (^2 = 0.2) and allowed the system to 
choose q which minimizes the free energy. The results 
are quite similar to the case when q = 0. We see that the 
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FIG. 11: (Color online) The phase diagram for 82 = 0, when 
the system can choose the value of q. Incommensurate SDW 
order appears below T* — 0.56T a and leaves some parts of the 
FS ungapped, allowing for co-existing SC order, (a) T s /T c = 
3. The SDW+SC phase appears only in a small region near 
T c . At low T the system still undergoes a first order transition 
between commensurate SDW and SC states, (b) For larger 
T a /T c — 5 (weaker SC interaction) the co-existence region 
widens and extends down to T = 0. The q = SDW state 
has the lowest energy at T — for So < 0.195. 
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FIG. 12: (color online) Same as in Fig.[8ja) but now we allow 
the system to choose the value of q. The phase diagram from 
Fig. [8] is shown by dashed lines. A finite q emerges below 
a particular T and moves the co-existence region to larger 
do, together with the SDW-N transition. This broadens the 
co-existence region, and slightly changes the shape of T c (So) 
inside the magnetic dome. 



SDW and SC orders do co-exists in the parameter range 
which extends from the crossing point down to T = 0. 
The width of the co-existence region widens a bit when 
we allow the system to choose q, but qualitatively, the 
behavior in Figs. 151 and H121 is the same. Note, in our two- 
band model, the ellipticity of of the electron FS breaks 
the rotational symmetry and favors the direction of q 
along the ellipse's major axis, see Eq. (|3.22b '). 

To summarize, SDW and SC -1 phases do co-exist in a 
range of finite dopings, but the width of the co-existence 
region depends on the amount of ellipticity of the electron 
band and the ratio of T s /T c . At larger T s /T c the width 
of the co-existence region increases for fixed 62 , and there 



is optimal 82 at which the width is the largest. The fact 
that the system can lower the energy by making SDW 
order incommensurate also acts in favor of co-existence, 
but qualitatively the picture remains the same as in the 
case when q is set to be zero. 



B. 



Minimal co-existence with s ++ state 



We next look at the SC state with gaps of the same 
signs on two FSs. Such states seem unlikely for pnictides, 
because they require a negative sign of the inter/band pair 
hopping termJ^ Still, it would be interesting to investi- 
gate consequences of attractive SC interaction between 
electron and hole bands. 

The expressions for £± in this case is slightly more 
complicated and less illuminating than those for 
state, although quite similar, and so are the self- 
consistency equations, which we do not write here, but 
which are obtained from Eqs. (|2.41j) - (|2.43[) in a way com- 
pletely analogous to Eqs. (|4.ip . We first present the re- 
sults for ^2 = 0, Fig. [13] We found that co-existence 
region does not appear even if we allow SDW order to 
become incommensurate. There are commensurate and 
incommensurate SDW phases on the phase diagram, and 
SC ++ phase, but the transition between SC and SDW 
phases remains first order. In other words, the appear- 
ance of gapless excitations in the SDW phase due to in- 
commensuration at large So does not seem to favor a 
mixed superconducting and magnetic state, in sharp con- 
trast to the case of SC, where incommensuration 
induces co-existence, see Fig. ITTTb). 

For a non-zero 62 , there might appear a tiny region of 
co-existence at low temperatures. We illustrate this in 
Fig-HH where in panel (a) we plot the phase diagram for 
82 = 0.2 and set q = 0. (When the system is allowed to 
choose q, the results change minimally, in a way similar 
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FIG. 13: (Color online) The phase diagram for a conventional 
s ++ SC order parameter, at 82 = 0, and varying So- We allow 
the system to choose q. The SDW+SC state does not appear, 
even when SDW order becomes incommensurate. 
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FIG. 14: (Color online) (a) Same as in Fig. [13] but for fixed 
82 = 0.2 and 5 = 0. For most of the phase diagram, the 
behavior is the same as for 82 = 0, but there appears a very 
tiny range of SC + SDW phase at the lowest T. The transi- 
tion to purely SC state is first order at all T. (b) The zero- 
temperature phase diagram. SDW and SC s ++ states are 
separated by first order transition virtually everywhere ex- 
cept a small region at finite So and 82, where SDW+SC state 
emerges. 



to Fig. fT2]) . In panel (b) of this figure we show where the 
region of SDW+SC ++ exists for different So and 82- We 
see that the range of co-existence is very narrow, and we 
also found that the difference in free energies between a 
pure SDW state and SDW+SC state is very small due to 
small value of the SC order parameter. 

Observe also that the co-existence region in Fig. Q3] is 
to the left of the line Sq + 82 = too at which a Fermi sur- 
face appears in the SDW state (a dashed line in Fig.H4b). 
In other words, s ++ superconductivity does not emerge 
even when there is a Fermi surface in the SDW state. 
This shows once again that the presence or absence of the 
Fermi surface in the SDW state is not the primary reason 
for the presence or absence of the SDW+SC phase. The 
true reason is energetic - the SDW+SC state can either 
lower or increase the energy compared to pure state de- 
pending on whether SDW and SC orders attract or repel 
each other. The absence of the co-existence phase even in 
the range where SDW state has a Fermi surface is a clear 
indication that there is the "repulsion" between SDW 
and SC orders, if the SC order is s++, Eqs. flggQ - lggS) . 
The same conclusion was recently reached by Fernandes 
et d& 



V. SDW + SC, ANALYTICAL RESULTS 

We corroborate the numerical analysis in the preced- 
ing Section with the analytical analysis. We first present 
the results of Ginzburg-Landau (GL) description near the 
point where second-order SDW-N and SC-N transitions 
meet, then consider the phase diagram at T = 0, and 
finally combine the two sets of results and compare ana- 
lytical phase diagram with Fig. [7] 



A. Ginzburg-Landau analysis 

We begin with the GL analysis near the point where 
T S (S ,5 2 ) = T c . Near this point, both the SDW and SC 
order parameters are small and we can expand the free 
energy, Eq. (|2.43j) . to the fourth order in to and A and 
compare different phases. For simplicity, in this section 
we assume that the SDW order is commensurate. An 
extension to a finite q complicates the formulas but does 
not change the outcome. 

The expansion of the free energy, Eq. (|2.43|) in powers 
of m and A yields 



J- = «aA 2 + a, 



+ ylA 4 + £m 4 + 2CA 2 



(5.1) 



where T = AF(m, A)/(4iVp). Coefficients oa, ce m , 
and B in Eq. (|5.ip are identical for both s" 1 and s 
SC states: 



.4. 




and 



A = 



B = 



irT v ^ 1 




(5.4) 
(5.5) 



The difference between s + and s ++ SC orders appears 
only in the coefficient C. For state we have 



C (+-) 




2\2 



while for s ++ 

a 



(++) 




4 ^W^^, 



(5.6) 



(5.7) 



Note, that, although both C-cocfficicnts are positive, 
this does not preclude co-existence in Eq. (|5.1[) . and we 



find below that the sign of parameter x = AB — C is 
more important for co-existence. We will demonstrate 



that since C, 



(++) 



> C(_| ), x i s positive for a broader 

range of parameters in s" 1 state than that in s ++ state. 
In fact, x remains always negative in s ++ state. Below we 
will use the notion that x > corresponds to an effective 
attraction between the two orders. 

The free energy, Eq. (|5.1[) . has two local minima, cor- 
responding to pure states, when one of the order param- 
eters is identically equal to zero: 

1) a pure SC state, defined by to = and dF/dA = 0, 
has the free energy and SC order parameter 



•Fa 



"A 

4 A'' 



A 



«A _ 

2A ; 



(5.8) 
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2) a pure SDW state, defined by A = and dJ- '/dm = 0, 
has the free energy and SDW order parameter 



2B' 



(5.9) 



In addition, the free energy may also have either a 
saddle point or a global minimum when both A ^ and 
m^O. To sec this, we write the free energy Eq. (|5.1[) in 
equivalent form, 



T = 



C a2 



m 2 + ^A 2 



Ot A - ~j^ a ™ ) A 2 



r<2 

A- — ) A 4 
B 



which is now a sum of two independent parts for A 2 and 
A/ 2 = m 2 + (C/B)A 2 . For an extremum state, given by 
8aJ~ = d m F = 0, the stationary values of order param- 
eters, 



A 



aA-B — a m C 



C 



determine the free energy, 



m&A 



-BM 4 - AB C V . 
B 



(5.10) 
(5.11) 



When both coefficients in Eq. (|5.11[) arc positive, 

B>0 , x = AB-C 2 >0, (5.12) 



the mixed state, Eq. (|5.10p . corresponds to the minimum 
of the free energy, which is smaller than the minima for 
pure SC or SDW states: 



J~ m&cA — J~ m 



Ta~ 



1 {aAB-a m C) 2 
AB AB - C 2 

1 (a m A-a A C) 2 
AA AB-C 2 ' 



(5.13) 



Consequently, in the phase diagram, the pure SDW 
and SC states are separated by a SDW+SC phase, and 
the transitions into this intermediate state are second- 
order. However, if B > and x < 0, the mixed phase, 
Eq. (|5.10[) . corresponds to the saddle point of the free 
energy and is not thcrmodynamically stable phase. In 
this case, pure SDW and SC phases are separated by a 
first-order transition line. When B < 0, one needs to ex- 
pand further in m to determine the phase diagram. We 
will not discuss the case B < further within GL theory. 

We apply Eq. (|5.12|) to the case 8^ = So + 82 cos 2<j> 
which we considered in the previous Sections. We remind 
that 82 = corresponds to co-circular FSs with different 
chemical potentials, while <5 = corresponds to FS ge- 
ometry in which k c F = hp, but the electron pocket is 
elliptical. 

At perfect nesting 5q = 82 = 0, and the system devel- 
ops an SDW order at T s > T c . Deviations from perfect 



nesting lead to two effects. First, as we already said, 
the magnitude of a m is reduced because SDW instability 
is suppressed when nesting becomes non-perfect. Super- 
conducting a a is not affected by 5^, and eventually wins 
over SDW. Second, coefficients B and C evolve with 
and, as a result, the sign of % = AB — C 2 depends on 
values of i5o and 82- 

The GL expansion is applicable only in the vicinity 
of points at which the temperatures of the SDW-N and 
SC-N transitions coincide T s (6^) = T c . This condition 
together with Eq. (j3.6[) establish the relation between Sq 
and 82 at which one needs to compute the parameters B 
and C. 



1. s + superconductivity 

To get an insight on how \ evolves with £r , we 
first assume that T s /T c is only slightly larger than one 
(T s /T c = l + 8t), in which case T s (8^) = T c at small S 
and 82, and we can expand A, B, and C in powers of 8q 
and 82- Specifically, we have from Eq. (|3.6[) 



0.663 



*o 2 



-A 



(5.14) 



where £(3) is a Riemann Zcta function. Collecting terms 
up to the fourth order in the expansion, we obtain 

where 

(§<sW) ■ 

The sums are expressed in terms of the Riemann-Zeta 
function C(3), C(5), and C(7) and give s\ rts 5.261 and s 2 ~ 
9.082. Substituting 8y. = 80+S2 cos 2<j> and averaging over 
momentum direction d> on the FSs. we obtain 



1 



X 



327r 8 ri* 



-3.820(5^ + 6.703^5| - 0.297<5|) . (5.17) 



We see that for 80 = 82 = 0, x = 0, i.e., for a perfect 
nesting the system cannot distinguish between first order 
transition and SDW+SC phase. This result, first noticed 
in Ref. [29T implies that the phase diagram is quite sen- 
sitive to the interplay between So and 82- We see from 
(|5.17[) that in the two limits when either 82 = or So = 0, 
X < 0, i.e., the transition is first order. This agrees with 
the numerical analysis in the previous Section. We em- 
phasize that in both limits, a small SDW order, which we 
consider here, still preserves low-energy fermionic states 
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FIG. 15: (Color online) The three-dimensional plot of the 
SDW-SC crossing surface, T s (So,S2) = T c . At each point on 
the surface we show the sign of B and x = AB — C 2 . In 
the region where x > 0, the transition between SDW and 
SC states occurs via the co-existence region. For x < 0, pure 
SDW and SC states are separated by the first-order transition. 
When B < 0, the SDW-N conversion is of the first order and 
the present GL analysis is invalid. 



near the modified FSs. Fermions near these FSs do have 
a possibility to pair into s H state. However, SDW+SC 
state turns out to be energetically unfavorable. We par- 
ticularly emphasize that the ellipticity of electron disper- 
sion is not sufficient for the appearance of the SDW+SC 
phase near T c ~ T s . 

When both Sq + 1 and 82 7^ 0, there is a broad range 



0.765 < -± < 4.689, 
8q 



(5.18) 



where x > and the transition from a pure SDW phase 
to pure a SC phase occurs via an intermediate phase 
where the two orders co-exist. This also agrees with the 
numerical analysis (see Fig. [5] and Fig. fTU|) . 

Eq. (|5.18[) has to be combined with the equation for 
T s {Sq,S2) = T c , and the boundaries in Eq. (|5. 18|) set 
the critical values of 62 and Sq as functions of T s /T c . 
Combining Eqs. f|5 . 1 8[) and (|5.14[) , we obtain that co- 
existence occurs for 



(5.19) 



0.826TO V<ft < 6 2 < 1.663m V(5t 



To verify that this result holds at larger values of So 
and 82, we computed x without expanding in 8^ . We 
plot the resulting phase diagram in Fig. [15] The result is 
qualitatively the same as Eq. (|5.17|) : for Sq — or 8 2 = 0, 
X < and the transition between SDW and SC states is 
of first order, while when both <5o and 62 are non-zero, 
there exists a region where x > and the transition from 
SDW to SC state occurs via an intermediate SDW+SC 
phase. 
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FIG. 16: (Color online) The case of s ++ SC. Panel (a): the 
behavior of x(8o, 5 2 ) = AB - C 2 for different T s /T c for three 
cases: x(8,0), x(0, <5), and x(8,8). For each case, 5 is chosen 
to satisfy the condition T 3 (So,S2) = T c for a given T s /T c . We 
see that in all three cases, x < no matter what the ration 
T 3 /T c is. Panel B - the coefficient _B(<5o,<52) along the line 
T 3 (5o,S2) — T c . GL analysis is only valid when B > 0. 



2. 



superconductivity 



We performed the same calculations for a conventional, 
sign-preserving s-wave superconductivity. The key dif- 
ference with the s" 1 case is that now x = AB — C 2 is 
non-zero already when So = 62 = 0. Substituting A and 
B from (I5~4l) and (l5~5l) and C from (IBTfl) we obtain 



7C(3) 
1287r 4 T4 



< 



(5.20) 



The implication is that, for small Sq and S2, X remains 
negative and the transition between SDW and SC states 
is first-order. This result was first obtained by Fernandcs 
et al. in Ref. [29|. These authors also argued, based on 
their numerical analysis of the free energy, that there 
is no SDW+SC phase for s ++ gap even when Sq = S 2 
are not small. We analyzed the sign of x f° r larger i5o 
and 82 using our analytical formulas and confirmed their 
result. In Fig.[T|)]we show the behavior of x(^o, 82) at the 
transition point T s ((5o, £2) = T c for three representative 
cases: x(^o>0), x(0i^2), and x{5o,8 2 = Sq). In all cases, 
when B > (e.g., our GL analysis is valid) x(^0j^2) 
remains negative. 

We caution, however, that the absence of co-existence 
between s ++ SC and SDW states within GL model does 
not imply that the two states are always separated by 
first-order transition. GL analysis is only valid near 
T s (8o,82) = T c , when both orders are weak. The situ- 
ation at lower T has to be analyzed without expanding 
in to and A. And, indeed, we did find a small co-existence 
region T = 0, see Fig. Q3] 
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B. Zero-temperature limit 

We consider only the case of SC and the limit when 
relevant So and 82 arc small, i.e., when T s /T c = 1 + St 
and St <C 1- We compare energies for pure SDW and SC 
state and for the co-existence state and find the region 
where the co-existence state is energetically favorable. 

For this, we first verified that, at small Sq and 82, the 
values of SDW and SC order parameters at T = remain 
the same as at 62 = Sq = 0, i.e., to = mo = 0.28073 x 
(2nT 8 ) and A = A = m (T c /T s ). These values only 
change at large enough <5o and 82, e.g., to changes when 
So + S 2 > m . 

The free energies of pure SDW and SC states for 
to, A > Sq + S2 can be straightforwardly evaluated at 
T = by replacing the frequency sums in (|2.43|) by inte- 
grals. We obtain 

_ . to 2 S 2 Sn m 2 , to . „ , 

Hm) = -- + f + f + -m_ (5.21) 

A 2 A 2 / A \ 
^(A) . -_ + _,„(_) (5.22) 

These free energies have minima at to = too and A = Ao , 
respectively. At the minima, 

m 2 S 2 S 2 

^(Ao) = -X = -(!) 2 X- (5 ' 23) 

Observe that .F(Ao) < J 7 (mo) when T c = T s . This is 
the consequence of the fact that SDW magnetism is de- 
stroyed by doping and ellipticity, while superconductivity 
is unaffected. 

Comparing J- (mo) and .F(Ao), we find that the first 
order transition between pure SDW and SC states occurs 
at 

m 2 St = S 2 + S 2 2 /2. (5.24) 

If there is no intermediat e co-existenc e phase, the SDW 
state is stable for Sq < Stm^ — <5 2 /2, while SC state is 
stable for larger values of Sq. 

We next determine when the intermediate state ap- 
pears at T = 0. For this we expand the free energy near 
the SDW and SC states in powers of A and to, respec- 
tively. We then obtain, near the SDW state, 

F(m, A) = J-(mo) + a A A 2 + fe A A 4 , (5.25) 

and near the SC state 

T(m, A) = J"(A ) + a m m 2 + b rn m 4 . (5.26) 

We verified that 6a and b m are positive, while a m and 
cja can be of either sign. The key issue is what are the 
signs of a m and <xa at the point where J- (mo) = J-~(Ao). 



We found that, to leading order in St, a m = aa = a at 
this point, and a is given by 

«=f (1-8. + 7, 2 ), z=^ (5.27) 

Note that to obtain a we had to expand to order St 2 . By 
virtue of Eq. (|5.24j) . Sq = m ^/St^l — z, i.e. we have to 
consider z < 1. 

When a is positive, both pure states are stable, and 
there is a first-order transition between them. When a < 
0, the pure SDW and SC states are already unstable at 
the point where J-(mo) = J-"(Ao), what implies that when 
we vary (5o at a fixed ($2, there is a range of Sq around 
Sq = mo's/Sts/l — z in which the co-existence state has 
a lower energy than the pure states. From (|3.1|) we see 
that a > when z < 1/7, while a < for 1/7 < z < 1. 
In terms of 62 , this implies that the transition at T = is 
first order between pure states when S2 < 0.535m VSi, 
while at larger <5 2 , pure SDW and SC phases are separated 
along Sq line by the region of the co-existence phase. The 
width of the co-existence phase initially increases as 82 
increases, but then begins to shrink and vanishes when 
S2 approaches S2 — lA\AmQ\/~8t (z approaches one from 
below). At this point, the co-existence region shrinks to a 
point Sq = 0. At larger S2, the SC state has lower energy 
than the SDW state for all values of Sq 

If we keep S fixed but vary <5 2 , the co-existence range 
appears at Sq = +0 (z = 1) and exists up to Sq = 
0.926to o v / ^ (z = 1/7). At larger S (z < 1/7), there is a 
first order transition between pure SDW and SC states. 

C. The phase diagram 

We now combine the results of GL analysis near the 
crossing point and at T = into the phase diagrams. 
For definiteness, we set St = T s /T c — 1 to be small and 
consider the set of phase diagrams in variables T and Sq 
for different fixed 82- The results of this subsection has 
to be compared with the phase diagrams presented in 
panels (al)-(a3) in Fig. [3 see also Fig. [T0l 

From the analysis in the preceding two subsections, 
we found five critical values of 82 ■ two are obtained from 
the GL analysis of the range of the co-existence phase, 
and are given by (|5.19p . two are critical values at which 
the co-existence phase first appears and then disappears 
at T = 0, and the last one is the maximum value of 82 
at which T S (S = 0,<5 2 ) = T c . From ([534]) this value 
is 82 = 1.739movSt. Arranging these five values from 
the smallest to the largest, we obtain the following set of 
phase diagrams at small St: 

(a) For 82 < 0.535too\/^, there is no intermediate 
phase, and pure SDW and SC transitions are sepa- 
rated by a line of a fist-order transition. The line is 
tilted towards smaller So at smaller T: it originates 
at So = m a VSi(1.508 - 0.5(5|/(m^t)) 1/2 at T = T c 
and ends up at So = mo\/~8i.(l - Q.hSlKm^St)) 1 / 2 . 
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FIG. 17: (color online) A partially-gapped SDW state, where 
only a fraction N„dw of the electron and hole FSs is nested. 
On the remaining parts the dispersions and £ c are very 
different. The SDW state appearing below T s gaps excitations 
only in the shaded/boxed areas (see Fig. Eld)), while on the 
rest of the FSs the dispersions are close to the original and 
£ c . The SC state below below T c does not compete with SDW 
in non-nested regions, but competes with SDW state in the 
nested (boxed) regions. 



(f) $2 > 1.739moV^t) T s (5q,52) becomes smaller than 
T c for all So, and the system only develops a SC 
order. 

This behavior is also totally consistent with Fig. [10) all 
different phase diagrams are reproduced if we take hori- 
zontal cuts at different 62- We see therefore that numeri- 
cal and analytical analysis is in full agreement with each 
other. 

The only result of numerical studies not reproduced in 
small St analytical expansion is the existence of a range of 
62 where the transition between the SDW phase and the 
co-existence phase is second order, while the transition 
between the SC phase and the co-existence phase is first 
order, see Fig. [TO] To reproduce this effect in analytical 
treatment, we would have to expand to the next order 
in St. Note in this regard that it is evident from Fig. 
ITUl that the width of the range where one transition is 
first order and another is second order shrinks as T s /T c 
decreases. 



(b) For 0.535m \/8t < S 2 < 0.826m VSt <, the inter- 
mediate phase appears near T = and extends to 
some T < T c . At larger T, the transition remains 
first order. This behavior is consistent with the 
panel (al) in Fig. [7J 

(c) For 0.826m o y/Si < S 2 < lAUmoVSt. the interme- 
diate phase occupies the whole region T < T c . This 
behavior is consistent with the panel (a2) in Fig. [7J 

(d) For lAUmoVSt < 5 2 < 1.663m Vsi, SC state 
wins at T = for all Sq. There is phase transi- 
tion at finite T. The transition is first order be- 
tween pure SDW and SC state at smaller T, but 
the co-existence phase still survives near T c . This 
behavior is consistent with the panel (a3) in Fig. [7J 

(e) For 1.663m V^ < S 2 < 1.739m Q VSt, the co- 
existence phase near T c disappears, and the transi- 
tion becomes first-order along the whole line sepa- 
rating SDW and SC states. 



VI. PARTIAL SDW STATE 

In previous Sections we considered the situation when 
the splitting between hole and electron FSs is small. We 
now consider how the phase diagram is modified if in 
some fc-regions hole and electron FSs are quite apart from 
each other (after we shift the hole FS by Qo). Such re- 
gions are far from nesting and we make a simple assump- 
tion that they are not affected by SDW. We then split 
the FS into nested parts where commensurate SDW state 
exists and a SC order can exist as well, and non-nested 
parts, where only SC order is possible. We present this 
schematically in Fig. 1171 The nested parts lie in some in- 
tervals of angles </> with total circumference A^>, and have 
weight N sdw < N total = 1 (A0/2tt = N sdw / N to tal) ■ 

The free energy and the self-consistency equations then 
can be written as sums of the two contributions. The first 
sum is over the FS part that has only SC order parameter, 
and in the second sum we integrate over part of the FS 
with both orders.— 
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FIG. 18: (color online) (a) The phase diagram for s H su- 
perconductivity, and SDW order parameter existing only in 
boxed regions of the FS in Fig. 1171 with the relative width 
N s dw We manually set the doping dependence of N s d w to 
be N s dw = 0.8 — 0.4(5o/0.3) and neglected the effect of this 
variation of N 3 dm on T s . Observe that SC and SDW orders 
co-exist in a wide range of So and T. (b,c) the order parame- 
ters, A and m, and the free energy, F, as functions of So at a 
constant temperature, T/T c — 0.2 (b), and as functions of T 
at a constant So = 0.13 (c). As a function of So, SDW order 
parameter starts decreasing when SC order appears, and then 
jumps to zero and the system becomes a pure SC. 



The self-consistency equations (|6.2|) and (|6.3j) are 
obtained by minimization of the functional AF, 
d(AF)/dA = and dF/dm = 0, and these expressions 
reduce to previous formulas (|2.41[) - ()2.43[) for N sdw = 1. 

We find that the results are very similar to what we 
found within the approximation of a small FS splitting. 
The typical picture is shown in Fig. [18] 

The only differences from Fig. [8] in this case are the co- 
existence of SC and SDW states already at zero doping 
Sq = 0, and weak first order transition to purely SC state. 
We also analyzed s ++ SC order and again found a much 
weaker tendency for co-existence, similar to Fig. 1141 



VII. CONCLUSIONS 

To conclude, we presented a general theoretical de- 
scription of the interplay between itinerant SDW and 
SC orders in two-band metals. Within the mean-field 
approach we derived coupled self-consistency equations 
for the order parameters and the expression for the free 
energy, which is necessary to determine the stability of 
different phases. 

We considered the FS geometry with one hole and one 
electron bands of different shapes (a simplified FS geome- 
try for Fe-pnictides) and investigated the phase diagrams 
and the stability of the SDW+SC states for: (a) differ- 
ent gap structures of the SC state, Figs. [51 [JTJ vs. [JJ] (b) 



variations in the relative strength of SDW and SC inter- 
actions, Figs. [T0l QTJ (c) cllipticity of electron pockets, 
Figs. [5J QJJ an d (d) incommensuration of SDW order, 
Figs. [Til EH We considered the case when the transi- 
tion temperature to pure SDW state, T s , is higher than 
the critical temperature T c of a pure SC state. In the 
opposite case, T s < T c , the SC state develops first and 
suppresses SDW state. 

We found that the SC s ± state with extended s-wave 
symmetry has much stronger affinity with the SDW state 
than the traditional s ++ state. A co-existence region of 
,s ++ SC state with SDW is tiny and the co-existence is 
anyway very weak in terms of energy gain compared to 
the pure SDW state. The transition from the pure SDW 
state to the pure SC state is always first order, Fig. [JJ] 
For s 1 * 1 gap, there is a stronger inclination towards co- 
existence with SDW state due to effective "attraction" 
between the two orders. We found that, depending on 
the interplay between different effects (e.g., cllipticity 
and doping), the transition between SDW and SC or- 
ders is either first order or continuous, via the interme- 
diate SDW+SC phase, in which both order parameters 
are non-zero, Figs. [H[5]£2 

We note that the co-existence region gets larger with 
increased strength of the SDW interaction relative to its 
SC counterpart, described by the ratio T s /T c . Thus gen- 
erally we should see better co-existence between SDW 
and SC states, if T s is increasingly larger than T c , Fig.QJJ 

Our results are in a disagreement with a common be- 
lief that, because SDW and SC states compete for the 
Fermi surface, the SDW+SC state should emerge when a 
pure SDW state next to the boundary of the co-existence 
region still has a modified Fermi surface at T = 0, and 
should not emerge when fcrmionic excitations in the pure 
SDW phase are fully gapped at zero temperature. We 
found that the key reason for the existence of the mixed 
SDW+SC state is the "effective attraction" between the 
SDW and SC orders, while the presence or absence of the 
Fermi surface in the SDW state at T — matters less. 
Specifically, we found cases when SDW and SC orders 
do co-exist even when fcrmionic excitations in the pure 
SDW phase are fully gapped at T = 0, Fig. [TUTa). and 
we also found, for s ++ pairing, that there might be no 
co-existence down to T = even when the pure SDW 
phase has a Fermi surface, Fig. [141 

The phase diagrams for s H gap are quite consistent 
with the experimental findings in pnictides. For example, 
first order transition in Fig. [6] looks very similar to phase 
diagram of 1111 materials (La,Sm)OFeAs, where FSs are 
more cylindrical. The co-existence region in Fig. [5] corre- 
lates well with doped 122 materials based on BaFc2As2, 
where hole and electron FSs are less nested. And Fig. [iJJ 
shows that one can get both SDW+SC phase and first- 
order transitions for the same SC state and the same 
family of materials. Our key result is that the way the 
doping is introduced into the sample will determine the 
nature of the FS changes, and the path it will take in the 
(<5oi (^-plane: whether through a first order transition or 
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through a co-existence region. In other words, we argue 
that there is strong correlation between how exactly FSs 
evolve upon doping and whether or not SC and SDW 
states co-exist. 

The final remark. In the literature, there exists 
a notion of "homogeneous" and "inhomogeneous" co- 
existence of SC and SDW orders. The latter is a 
metastable state when the two orders exist in different 
spatial parts of the material. What we emphasize is that 
the other kind, "homogeneous" co-existence of SC and 
SDW orders in real space, is in fact "inhomogeneous" 
in momentum space: the SC and SDW orders dominate 
excitation gaps on different parts of the FS. 



Appendix A: Electron and Hole dispersion for small 
FS splitting 



In this Appendix, we discuss in detail the approxima- 
tion we used for the dispersions of holes and electrons 
for the case when the splitting between hole and electron 
FSs is small. 

SDW and SC orders mix c-fermions with momenta k 
and /-fcrmions with momenta k + q. The generic expres- 
sions for the two dispersions are 
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When the two FSs are circles of non-equal size, mt x = 
mfy = irif, we have (/z c ,m c ) rj (/x/,m/) w (fi,m), but 
m c =/= rrif and fi c =/= fj,f. The approximation we used in 
the text implies that 
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Introducing fkq and (5 kq defined in (|2.30[) . we obtain 
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We emphasize that, within this approximation, £k q and other depends on the angle along the FS. When mf x ^ 
^kq arc ^ w0 independent variables, one depends on the m fyi the derivation remains the same but 5 kq acquires 
deviation along the FS in the transverse direction, an- an additional term with the difference between m f x and 
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mfy (the cos 20 term). 
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